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One suggestion for determining the properties of QCD at finite temperatures and den- 
sities is to carry out lattice simulations with an imaginary chemical potential whereby 
no sign problem arises, and to convert the results to real physical observables only 
afterwards. We test the practical feasibility of such an approach for a particular class 
of physical observables, spatial correlation lengths in the quark-gluon plasma phase. 
Simulations with imaginary chemical potential followed by analytic continuation are 
compared with simulations with real chemical potential, which are possible by using a 
dimensionally reduced effective action for hot QCD (in practice we consider QCD with 
two massless quark flavours). We find that for imaginary chemical potential the system 
undergoes a phase transition at \fi/T\ ~ n/3, and thus observables are analytic only 
in a limited range. However, utilising this range, relevant information can be obtained 
for the real chemical potential case. 
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1. Introduction 



Given the applications to cosmology and heavy ion collision experiments, it is impor- 
tant to determine the properties of QCD at finite temperatures and baryon densities. 
For instance, one would like to know the locations of any phase transitions, and the 
properties of the quark-gluon plasma phase such as its free energy density, or pres- 
sure, as well as the spatial and temporal correlation lengths felt by various types of 
excitations in the system. 

Because the theory is strongly coupled, the only practical first principles method 
available for addressing these questions is lattice simulations. While there has been 
steady improvement in the accuracy of results at vanishing baryon density M, the case 
of a non- vanishing density is still largely open, despite much work Hi]-||13|. Indeed, 



introducing a non-vanishing density, or chemical potential, is difficult because it leads 
to a measure which is not positive definite (this is the so called sign problem), whereby 
standard Monte Carlo techniques fail. 

In this paper we focus on one of the suggestions for how a finite density system could 
eventually be addressed with practical lattice simulations. The idea is to first inspect 
an imaginary chemical potential, whereby the sign problem temporarily disappears, 
and then relate this to the case of a real chemical potential. Let us denote by \x the 
chemical potential for quark number Q, and by the chemical potential for baryon 
number B = Q/3: then \i = /ie/3. By fiR,fii 6 R we denote the real and imaginary 
parts of /i: 

H = H R + iH I . (1.1) 

It is easy to see, by going to momentum space, that physical observables are periodic 
in hi with the period 27rT. 

There are then two types of suggestions for how an imaginary jj, — could be 
utilised to obtain information on a system with a real chemical potential, \i = fiR- 
The first idea is related directly to the equation of state, and employs the canonical 
partition function at fixed quark number (|fj, || [|, [j| and references therein): 

Z(T, Q) = —J^ dfxi Z(T, n = m) e~^ Q / T , (1.2) 

where the grand canonical partition function, 

Z(T,n) = Tre- ( ^-^ )/T , (1.3) 

has been evaluated with an imaginary chemical potential. Here H and Q denote the 
Hamiltonian and quark number operators, respectively, while Q is a number. With an 
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imaginary chemical potential, Z(T, fi = ifij) or rather the ratio Z(T, [i = i/ii)/Z(T, 0), 
can be determined using standard lattice techniques. What remains is to perform the 
integral in Eq. ( |1.2| ). Of course, this gets more and more difficult in the thermodynamic 
limit Q — > oo, because oscillations reappear in the Fourier transform. In addition, a 
Legendre transform would be needed to go from Z(T, Q) to a system in an ensemble 
with a real chemical potential, Z(T,fi = [Mr). 

The second idea ([11] and references therein) is that, away from possible phase 
transition lines, the partition function and expectation values for various observables 
should be analytic in their arguments, in particular in fi/T. Thus, we may attempt 
a general power series ansatz for the functional behaviour in fi/T, determine a finite 
number of coefficients with an imaginary chemical potential, and finally analytically 
continue to real values. 

In this letter we study the latter suggestion within the framework of a dimensionally 
reduced effective theory. As we shall review in the next section, at temperatures 
sufficiently above the phase transition, the thermodynamics of QCD can be represented, 
with good practical accuracy, by a simple three-dimensional (3d) purely bosonic theory. 
This can also be done with a chemical potential, both real and imaginary ||14|| . We 
then use this theory to measure the longest static correlation lengths in the system 
for both cases. We find that, for small \(J*/T\, the observables are well described by 
a truncated power series with coefficients determined by fits. We then inspect how 
well the analytically continued series describes the real data. In principle the free 
energy density could be addressed with similar effective theory methods fl5fl , but this 
requires a number of high-order perturbative computations which are not available at 
the moment for fi ^ 0. 



2. Effective theory 
2.1. Action 

The effective theory emerging from hot QCD by dimensional reduction ||16|]-||22|], JT3 



is the SU(3)+adjoint Higgs model with the action 

S = Jd 3 x j^Tr 4 + Tr [A, A ] 2 + m 2 3 Tr A 2 + i l3 TrA 3 + A 3 (Tr^) 2 } , (2.1) 

where Fij = diAj — djAi + ig^Ai, Aj], Di = di+ ig^Ai, Fy, A i: and A are all traceless 
3x3 Hermitian matrices (Aq = A^T a , etc), and gf and A3 are the gauge and scalar 
coupling constants with mass dimension one, respectively. The physical properties of 
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the effective theory are determined by the three dimensionless ratios 



x= h nflfe = »8) I= T| (2 . 2) 

#3 #3 #3 

where /X3 is the MS dimensional regularization scale in 3d. For vanishing chemical 
potential 73 = and no term cubic in Aq appears. These ratios are via dimensional 
reduction functions of the temperature T / Ay^ and the chemical potential fi/A^g, as 
well as of the number Nf of massless quark flavours; for the case \x = 0, we refer to p2 



The inclusion of quark masses is also possible in principle, but in the numerical part 
of this work we assume Nf = 2 massless dynamical flavours, the other flavours being 
approximated as infinitely heavy. 

The mass parameter 777J, represented by y in Eq. fl2.2|) , turns out to be positive p2 



This guarantees that the 3d theory tends to live in its symmetric phase, Aq ~ 0, at 
least on the mean field level. We will return to this issue presently. 

Compared with the case \i = 0, the dominant changes in the action due to a small 



chemical potential are now 14 



(,3) 

where iV c = 3. Thus, one new operator is generated in the effective action, and one of 
the parameters which already existed, gets modified. 

For real chemical potential, \x = fiR, the effective action is thus complex, whereas for 
imaginary chemical potential, /i = ifij, it is real. 



2.2. Ranges of validity 

There are several requirements for the effective description in Eq. ( ]2.1|) to be reliable. 
They are all related to a sufficiently "weak coupling" , or effective expansion parameter, 
for a given T/Ajgg, /i/A^g, Nf. Let us briefly reiterate them here. 

First, the perturbative expansions for the effective parameters in Eq. (|2.2j) have to 
be well convergent. Inspecting the actual series up to next-to-leading order, it appears 
that this requirement is surprisingly well met even at temperatures not much above 
the critical one |£J . 

Second, the higher dimensional operators arising in the reduction step which are 
not included in the effective action in Eq. (|2.1|) , should only give small corrections. 
This condition is met if the dynamical mass scales described by the effective theory 
are smaller than the ones ~ 2nT that have been integrated out. In pure Yang-Mills 
theory, there is evidence that this can be sufficiently satisfied at temperatures as low as 
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T ~ 2Aj^g [Jr7|,[22]-|26[|,[r4"|. However, when fermions are included and a real chemical 
potential is switched on, some of the mass scales increase (see below), and the effective 
description will become less accurate. 

Third, the effective 3d theory represents the 4d theory reliably only when it lies in its 
symmetric phase |27j (i.e., A ~ 0). Indeed, for Nf = QCD has a so called Z(N)- 



symmetry p8| |29| , and this symmetry is not fully reproduced by the effective theory. 
The 4d Z(iV) symmetry is however spontaneously broken for Nf = in the deconfmed 
phase, and even explicitly broken for Nf > 0. In this case broken Z(iV) means that 
the Polyakov line is approximately unity, corresponding to A ~ and hence to a 
symmetric phase in terms of the gauge potential. Consequently the requirement to be 
in the symmetric phase of the 3d theory is easier to control for Nf > 0. Furthermore, 
the situation gets even better for real /j ^ 0. In the effective theory this can be 
seen, for instance, from the fact that the mass parameter y in Eq. (|2.3| ) grows, which 
makes it more difficult to depart from Aq ~ 0. Another stabilising factor is that the 
unphysical minima correspond to non-zero expectation values for TrAfj [p7l, and the 
imaginary term ~ zTr^Q in the action in Eq. (|2.1| ) disfavours such minima according 
to the standard argument [[3(J . 

On the other hand, an imaginary chemical potential fi = ifij, favours those Z(iV) 
broken minima where the Polyakov line has a non-trivial phase, and correspondingly the 
gauge potential is non-zero, A ^ 0. Utilising the perturbative effective potential J3~T| , 
we find that the lowest such minimum becomes degenerate with the symmetric one 
A ~ already at fii/T = tc/3, and increasing \ii further it eventually becomes lower 
than our minimum. Thus there is a (first order) "phase transition" ||. In the effective 
theory, this phase transition is triggered by the decrease of the mass parameter y in 
Eq. ( |2.3|) . Moreover, in this case the term ~ Tr A\ in the action in Eq. ( [2.1|) favours 
the effective theory remnant of one of the minima with Aq ^ 0, i.e., a non-trivial phase 
of the Polyakov line. 

This phase transition limits the applicability of the effective theory with imaginary 
chemical potential, since only the symmetric phase is a faithful representation of the 
4d theory ft7|] . (For one suggestion on how perhaps to circumvent this problem at 
least for Nf = 0, which we shall however not dwell on here, see ||32|| .) In a way, this 
problem is related to the fact that as one approaches Hi/T = ir, fermions start to obey 
Bose-Einstein statistics and become "light" infrared sensitive degrees of freedom (see 
also |31|), whereby it is no longer legitimate to integrate them out. 

In summary, the effective theory roughly loses its accuracy with a real chemical 
potential once even the longest correlation length is shorter than ~ l/(2nT), and with 
an imaginary chemical potential once \fi/T\ exceeds unity. Fortunately, this range of 
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validity contains the parameters that are phenomenologically most relevant. Indeed, 
heavy ion collision experiments at and above AGS and SPS energies can be estimated 
to correspond to ^b/T <4.0 [3~3"|, or a quark chemical potential \xjT < 1.3. 



2.3. Observables and their parametric behaviour 

As we have mentioned, the physical observables which we shall study are spatial cor- 
relation lengths: we consider operators living in the (x±, X2)-plane, and measure the 
correlation lengths in the ^-direction. 

In the presence of \x 7^ 0, there are only two different quantum number channels to 
be considered, distinguished by the two-dimensional parity P in the transverse plane. 
The lowest dimensional gauge invariant operators in the scalar (J = 0) channels are: 

J p = 0+ : Tr Al Tr i^, Tr A 3 , Tr A F? 2 , ... 

J p = 0" : TrFf 2 ,Tr^F 12 ,TrAoF 12 , ... (2.4) 

The corresponding 4d operators can be found in [R|J]. We shall measure whole cross 
correlation matrices between all (smeared) operators in these channels, but mostly focus 
on their lowest eigenstates, corresponding to the longest correlation lengths in the 4d 
finite temperature system. We denote the "energies" of these eigenstates, viz. inverses 
of correlation lengths, by m. We also examine the overlap of operators of different field 
contents onto the eigenstates. 

Since a change fi — > — fi can be compensated for by a field redefinition A — > —A in 
Eq. ( |2.1|) , all physical observables must be even under this operation. In the original 
4d theory the same statement follows from compensating /i — > —fi by a C (or CP) 
operation. Moreover, since there are no massless modes at /i = 0, we expect the 
masses to be analytic in fi away from phase transitions. For small values of fi/T, the 
inverse correlation lengths may thus be written as 

We have chosen to include 7rT in the denominators, because the chemical potential 
appears with this structure in the effective parameters, cf. Eq. ( |2.3| ). Of course, the 
radii of convergence of such expansions are not known a priori. 

Here we first check to what extent a truncated series of the type in Eq. ( |2.5| ) can 
accurately describe the data. In the range where this is possible, we determine the {q} 
with /1 = ifij, and check if the analytically continued result reproduces the independent 
measurements carried out with n = ji R . 
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Table 1: 

determines the lattice spacing (for the detailed relations employed here, see WM). 



The parameters used for |i / (cf. Eq. ffli.ty ). All correspond to T = 
--2. In addition, x = 0.0919, g\ = 2.92T, (3 = 21, volume = 30 3 ; where (3 



3. Simulations 



3.1. Simulation methods 

We simulate the theory at several yu/T. The values chosen, together with the cor- 
responding continuum parameters, are listed in Table |l[ Discretization and lattice- 



continuum relations |35[ are implemented as in [[14]]. As discussed there, finite volume 
and lattice spacing effects are expected to be smaller or at most of the same order as 



the statistical errors for the parameter values we employ. Compared with |L4j], we have 
increased the statistics and included many new values of fi/T, in order to carry out 
more precise fits. 

For real fi = /i/j, the action in Eq. ( pTT| ) with parameters as in Eq. ( [2.3| ) is complex, 
which precludes direct Monte Carlo simulations. We must thus carry out simulations 



using a reweighting technique, which has been explained in detail in ||14|| . There it was 
found that physically realistic lattice volumes may be simulated for chemical potentials 
up to fiji/T<4. For imaginary ji = z/ij, the action in Eq. (|2.1|) with parameters as in 
Eq. (|2.3j ) is real, and correspondingly we simulate the full action using a Metropolis 
update. 
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Figure 1: Inverse correlation lengths in the channel + , for real fi/T. "Scalar" states 
(Tr Aq etc) do depend on \ijT , while "gluonic" states (Tr F± 2 etc) are practically inde- 
pendent of it. For comparison, the horizontal band indicates the 3d pure glue result for 
Tr F± 2 W% i converted to our units via g\ = 2.92T. 



3.2. Results 

As a first result, let us note that, as has been the case in several related theories IBB" 



24], |14|, we again observe a dynamical decoupling of operators, such that operators 
involving scalars (Tr Aq, Tr AoFf 2 etc.) and purely gluonic operators (Tr F? 2 etc.) have 
a mutual overlap consistent with zero. The correlation matrix thus assumes an ap- 
proximately block diagonal form. We find that the gluonic states remain extremely 
insensitive to fi/T, and agree well with the masses found in d = 3 pure gauge the- 
ory pTfl . This situation is illustrated in Fig. |I[ 

The scalar states, on the other hand, show a marked dependence on fi/T, with their 
masses increasing for real fi and decreasing for imaginary values. For both small real 
and small imaginary fi/T, the ground state in each channel is scalar in nature, and we 
plot these states in Fig. §. 

Because of the different qualitative behaviours of 3d gluonic and scalar states, we 
may expect to observe a change in the nature of the ground state excitation at some /ir. 
Indeed, Fig. [I] suggests a level crossing at fi/T ~ 4.0. This would mean that the longest 
correlation length in the thermal system does not get arbitrarily short with increasing 
density, but rather stays at a constant level. Note that the value oim/T at this crossing 
is already so large that the effective theory may be inaccurate quantitatively, and in 
fact in the full 4d theory the flattening off could take place much earlier. However, the 
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Figure 2: Left: Inverses of longest correlation lengths in the channel + , for real and 
imaginary /i. Right: the same for 0" . 
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Table 2: Fitting the lowest masses in the channel + from jj = up to (i = /x max . 
The numbers in parentheses indicate the error of the last digit shown, the coefficients 
refer to Eq. \2. b\ ), the sub and superscripts r,i denote real or imaginary fi, and Q is 
the quality of the fit. 



qualitative effect should be the same. 

Next, let us discuss the applicability of the power series ansatz in Eq. (|2~5|) . To this 
end we perform fits over a range = 0...;U max to the inverses of the longest correlation 
lengths, both for real and imaginary fi. For imaginary /i, we can follow the "analytically 
continued" metastable branch as long as tunnelling into an unphysical minimum does 
not become a problem, which in practice means fii/T< 1.5. 
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Table 3: Fitting the lowest masses in the channel . The notation is as in Table ^. 

The results are shown for the + channel in Table 0, and for the 0~ channel in Table [| 
Examining these fits we see that in all cases we have good fits, as demonstrated by the 
low x 2 /dof and good Q values. In the case of real \i we find stable and well constrained 
values for the coefficients as we increase the size of the fitting range. For imaginary 
/j,, due to the breakdown of the effective theory at large values of fii/T, we have fewer 
significant data points, and consequently the coefficient of the quartic term is much 
less constrained. 

As our main result, we can now state that we observe good evidence for analytic 
continuation in the first non-trivial term, with cf consistent with c[ in the + channel 
and similarly for the 0~ states. Unfortunately, the data is not accurate enough to 
make a similar statement for cf, c\. Extremely precise measurements would be needed, 
because the range in \\xjT\ available to imaginary chemical potential simulations is 
very limited. On the other hand, from the phenomenological point of view the first 
non-trivial coefficient is sufficient, since the series expansion turns out to be in powers 
of fi 2 /(irT) 2 , which is small in the most important practical applications. Thus, for 
phenomenological purposes, it does not seem necessary to invest an extra amount of 
effort on a more precise determination of the masses in the imaginary /i case. 

4. Conclusions 

In this letter, we have studied the question as to what extent imaginary chemical 
potential simulations could be useful for determining the properties of the quark gluon 
plasma phase at high temperatures and finite densities. The physical observables we 
have measured are static bosonic correlation lengths, but the pattern should be very 
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similar for the free energy density, as long as T > T c . 

The method we have used is based on a dimensionally reduced effective field theory. 
This way we can address both a system with a real and an imaginary chemical potential, 
as long as their absolute values are relatively small compared with the temperature. 
For larger absolute values of fi/T = ifij/T, there is a (first order) phase transition, and 
the effective description breaks down. 

Despite the fact that we are only working in the quark-gluon plasma phase, we find 
an interesting structure in the longest correlation length, which decreases first but 
becomes constant beyond some real value of fi/T, which we estimate to be < 4.0. 

Furthermore, in the region where the effective theory is applicable, we find that direct 
analytic continuation does seem to provide a working tool for determining correlation 
lengths. For phenomenological applications, only the first two coefficients in the power 
series are needed, since we find the expansion parameter to be < fi 2 / (nT) 2 , which 
is small in heavy ion collision experiments. This is good, since determining more 
coefficients with imaginary chemical potential would require very precise simulations. 

We are thus encouraged to believe that in 4d simulations analytic continuation of 
imaginary chemical potential results would give physically relevant results if a good 
ansatz for the /^-dependence is available, and would allow to go closer to T c determining, 
e.g., the free energy density and the spatial correlation lengths there. Furthermore, it 
appears that at high enough temperatures it is even sufficient to determine only the 
coefficients of the first terms depending on /i, which amounts simply to susceptibility 
measurements at // = (see, e.g., |38], |39| and references therein). 
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